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Abstract 

The two theories for magnetic reconnection, one of Sweet and Parker, and the other 
of Petschek, are reconciled by exhibiting an extra condition in that of Petschek which 
reduces his theory to that of Sweet and Parker, provided that the resistivity is constant 
in space. On the other hand, if the resistivity is enhanced by instabihties, then the 
reconnection rate of both theories is increased substantially, but Petschek's rate can 
be faster. A different formula from the usual one is presented for enhanced Petschek 
reconnection. 

1 Introduction 

As is well known, the process of magnetic reconnection is important in many space 
and astrophysical contexts. The initial problem that first inspired research into the 
subject was the solar flare phenomenon, in which it appeared that energy was first 
slowly built up and stored in the magnetic field, and then suddenly released into ther- 
mal and kinetic energy. The first solution of the problem was given independently by 
Sweet (1958) and Parker (1957), who approximated the problem as a two dimensional 
incompressible MHD problem. They showed that the problem was essentially a bound- 
ary layer problem, and they estimated the rate of reconnection from a boundary layer 
analysis. This boundary layer analysis led to release of magnetic energy over a period 
of time several orders of magnitude longer than the observed energy release time in 
solar flares. A probable explanation of this discrepancy could be the fact that their 
estimates of the reconnection rate are based on normal (Spitzer) resistivity, while in 
the actual solar flare the resistivity could be greatly enhanced, leading to a much faster 
energy release. 

On the other hand, at the time when Sweet and Parker developed their theories, 
the possibility of enhanced resistivity was not appreciated, and other means of increas- 
ing the reconnection rates were sought. Petschek (1963) pointed out that, since the 
magnetic reconnection was a topological process, the field lines need not reconnect re- 
sistively along the entire length of the boundary layer, but could merge over a shorter 
length L'. 
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For this to happen, the rest of the boundary layer region should consist of slow 
shocks that could accelerate the matter that did not pass through the diffusive region. 
He found that the resulting reconnection rate was increased by the factor ^JL/L' , and 
by choosing L' small enough, a very rapid reconnection could be achieved. 

Ever since the two theories, Swcct-Parkcr's and Pctschck's, were published, there 
has been a controversy over which one was the correct one to apply. The controversy 
seemed to be settled, by the rather complete numerical simulation of Biskamp (1986), 
to be in favor of the Sweet-Parker result. Since then, a number of numerical simulations 
have confirmed this. 

Since both theories seemed rather well founded, it is a question of how either of them 
could be incorrect. In this note, I will show that the Petschek theory, as he proposed, 
it is indeed not correct, at least in the context of MHD with constant resistivity. In the 
development of his theory, Petschek left out one condition that also must be satisfied. 
The satisfaction of this condition leads to a uniqiic determination of the length L' in 
his theory, and indeed, if the resistivity is constant in space, it is the case that V is 
equal to L. This reduces Petschek's enhancement factor \fLjTJ to unity and Petschek's 
reconnection rate to Sweet-Parker's rate for constant resistivity. 

On the other hand, if one considers the possibility of enhanced resistivity, two things 
happen. (1) The Sweet-Parker reconnection rate becomes much faster, for the solar 
flare case, and (2) because such enhanced resistivity is very sensitive to current density, 
it can be space dependent also. This leads to L', being smaller than L, and to a even 
faster, Petschek like, reconnection rate. 

2 The Boundary Layer 

Reconnection is as much a global phenomena as a local one. For example, consider the 
magnetic reconnection of two cylinders with opposite poloidal flux. (Figure 1.) Let us 
also assume that the velocities induced by magnetic reconnection are slow compared to 
the Alfven speed everywhere, except in the reconnection and separatrix layers. Then, 
everywhere else, we have 

j X B = Vp (1) 

Also, since the layers are thin, we have the jump in p + B'^ /%tt zero across these layers. 
This means that, if the amount of reconnected and unreconnected flux is given, and if 
the rotational transform and pressure are known on each magnetic surface as functions 
of the poloidal magnetic flux, then there is a unique equilibrium. But as magnetic 
reconnection proceeds, one can keep track of the pressure and the rotational transform 
in both the regions of reconnected and unreconnected flux. The two regions change 
geometrically and physically as reconnection proceeds. Thus, the plasma first moves 
from the unreconnected region into the reconnection layer, where it is heated, then it 
flows into the separatrix region, and, finally, as the magnetic configuration changes, 
into the reconnected region. (However, this does not change the uniqueness of the 
equilibrium at each stage of reconnection.) 

Because of this uniqueness, the length of the reconnection layer L is totally deter- 
mined at each stage as well as the horizontal field just outside of the layer. 

Appreciating this fact, all three authors took the length of the layer as well as the 
variation of the bounding field as given. They assumed that the region into which 
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Reconnection 
layer 

A(P+B2/87i)=0 



Figure 1: Slow merging of two cylinders. The situation is in a slowly evolving quasistatic 
equilibrium everywhere, except in the layers across which p + B^/Stt is continuous. R is the 
region of reconnected flux, and U is the region of unreconnected flux. 

the plasma flowed (the separatrix region) was at the same pressure as the upstream 
ambient pressure. This was the boundary layer problem to be solved. 

3 The Sweet-Parker and Petschek Theories for 
a Constant Resistivity 

Let us now examine the two theories. First consider that of Sweet and Parker. The 
reconnection layer is sketched in figure 2. It is easily shown that the flow out of 
the layer is at the Alfven speed, Vx = Va = Bq/^/Attp. The incoming flow of matter 
—LVy = LVr must balance the outgoing flow VaS, where Vr, the reconnection velocity, 
is the incoming velocity outside of the layer, where the plasma is tied to the field lines. 
6 is the half thickness of the layer. Thus, 

VrL = Va5. (2) 

On the other hand, by Ohm's law the field diffuses up stream with a velocity r]c/An6 
with respect to the incoming plasma, with velocity —Vr. Thus, in steady state, 

Vr = r]c/An5. (3) 




Separatrix 
A(p+B2/87i)=0 
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Figure 2: The Sweet-Parker layer. In a steady state the magnetic diffusion velocity rjc/AirS 
balances the incoming reconnection velocity Vr, and the inffowing mass ^VrL balances the 
outgoing mass 45 V^. 



Prom these equations we obtain 



and 
where 

r/c/47r 

This is the Sweet-Parker result. 

The Pctschek theory is indicated in figure 3. In this model the diffusive region, in 
which the merging actually takes place is of a much shorter, length L', than L. The 
remaining length of the boundary is occupied by slow shocks. In the diffusive layer 
the behavior is similar to the Sweet-Parker layer, the main difference being that the 
acceleration of the velocity up to the Alfven speed along the layer, is accomplished 
by magnetic tension associated with a transverse field component By. (In the Sweet- 
Parker theory this acceleration is produced mainly by a pressure gradient.) Outside 
of the Petschek diffusive layer the acceleration up to Va is accomplished almost in- 
stantaneously by the slow shocks. The Sweet-Parker model for their diffusive layer is 
replaced by the identical conditions for the Petschek model, but with L replaced by 
L', leading to the Petschek reconnection velocity, 




a factor of ^JL/L' faster than the Sweet-Parker reconnection velocity. The shocks in 
the outer L — L' region reduce the upstream to zero, and accelerate the plasma 
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Bx=Bo(1-x^L2) 




Figure 3: The Petschek model. The diffusion region of length L' < L is the same as figure 
2, with L taken equal to L'. Outside the diffusive region shocks accelerate the plasma and 
reduce to zero. The shocks' y velocity By/ y/Aiip balances the incoming Vr velocity in a 
steady state. The By field is vital to this model. 



crossing them to Va in the x direction to match the plasma flowing out of the diffusive 
region with the same Alfven x velocity. 

The shocks propagate in the y direction upstream into the plasma with velocity 
By/y/iFp. Since the plasma is flowing with the reconnection velocity Vy = —Va we 
have in steady state, 

-Vr, (8) 



which determines the magnitude of By the transverse field component. This y compo- 
nent of the field increases linearly along x in the diffusive region from to this value, 
and it turns out that the tension produced by the jzBy force is just enough to accelerate 
the plasma in the layer up to the Alfven speed. 

These results are all given in Petschek's paper and present a nearly complete, but 
qualitative, physical picture for magnetic reconnection, encompassing the possibility 
of a diffusive layer much shorter than L. Further, in his theory, L' appears to be a 
free parameter. Petschek then chooses L' as short as possible to get the maximum 
reconnection velocity. He determined this minimum to be the lower limit so that the 
current in the shocks did not seriously perturb the incoming magnetic field B^ = Bq. 
This limit was roughly 

L'>|(ln5)2, (9) 
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so that substituting in equation (7) he got a very fast Umit on the reconnection velocity 



Vn < ^. (10) 

This latter hmiting velocity has been generally quoted in the literature as the so-called 
Petschek reconnection velocity, and it is this velocity that has been compared with the 
Sweet-Parker reconnection velocity equation (4) in the controversy between the two 
theories. The disagreement should more appropriately be between the Sweet-Parker 
velocity and equation (7). 

Now, it turns out that L' is not a free parameter in Pctschck's theory. There is 
an additional condition associated with By that Petschek did not include . The all 
important By field, which is needed to support the shocks, is embedded in the rapidly 
moving plasma and is swept down stream at the Alfven velocity, Va- This field which 
is being swept away so rapidly must be regenerated at the same rate to preserve a 
steady state. 

This regeneration occurs during the merging process. The external field is nonuni- 
form, being strongest near a; = 0, so the lines of force at the center of the diffusive 
layer will move into the diffusive region fastest there. However, the nonuniformity of 
the external field is on the scale of the length of the layer L because of the global con- 
ditions, so as L' gets smaller, the nonuniformity over the shorter length is also smaller, 
and the regeneration process is weaker. A balance between the nonuniform merging 
process that creates the By field, and the down sweeping which destroys it, must be 
reached and this leads to a relation between L' and By. Thus, combining this relation 
with equations (7) and (8), determines L' uniquely. 

Let us estimate this balance qualitatively. The equation for By will be shown to be 

dBy_VnL'^ j.Va_ 

The second term on the right is the down sweeping term that destroys By . Its form is 
obvious. 

The first term represents nonuniform merging, and its form can be derived as 
follows: (see Figure 4). 

The external field depends on x as 

B^ = Bo{l-xyL^). (12) 

We assume that each fresh line that is merging enters the layer with velocity, 
proportional to Bx{x), so that the line enters faster at x = than at a; = L'. Thus, 
after entering the layer it will turn at the rate 

d9 _ V^{0)-V^{L') _ rjc , , 

di- L' /^^J' 

or 

-/2 



de rjc L'^ 1 



dt 4tt5 L2 L 



(14) 



The turning of a line of strength Bq at the rate d9/dt produces a component By at the 
rate {d6/dt)Bo = Vr{L''^ / LP')Bq/ L\ which gives the first term in equation(ll). 
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dBy/dt= (de/dt)Bo 



Figure 4: Regeneration of the By field by nonuniform merging. This rotates the field 
into the y direction. 



Now, setting dBy/dt = 0, for a steady state, gives from equation(9), 

By = ^^{L^/L^)B,. 



Prom equation (8) we get 

Vr = ^{L'-'/L') 

or. 



Va 



L' = L 



(15) 

(16) 
(17) 



Thus, L' is no smaller than L, and the Petschek rate equation (7) reduces to the 
Sweet-Parker rate equation (4). 

I believe that this is the reason that the numerical simulations always yield the 
Sweet-Parker rate, rather than the faster rate implied by Petschek's formula. 

A more formal derivation of equation (11) is given in the appendix. 



4 Anomalous Resistivity in the Sweet-Parker 
Model 

In the absence of a Bz component (that is no guide field), there is a strong instability, 
the lower hybrid instability, that should be excited, (Davidson 1975) This is the case if 
the current density in the layer is large enough that the difference in the electron and 
ion bulk velocities V- and Vj^ is greater than the ion acoustic speed. That is, if the 
drift velocity Vd = vj^ — v _ satisfies Vd > Vi with Vi the ion thermal velocity. Note that 
j = ne{e/c)vd. 
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There arc three wen-documented instances where magnetic reconnection is defi- 
nitely taking place, i.e. the solar flare, the magnetosphere-solar wind interface, and 
the magnetotail. If one examines these three cases, and applies the Sweet-Parker model 
to them, one finds in all three cases that the drift velocity, v^, is much larger than Vi. 
One can express this as follows: There is a critical current, jc, and critical layer thick- 
ness 5c = Bq/AttJc such that if 6 < 6c and therefore j is greater than the critical j, 
then the lower hybrid mode should be excited. This lower hybrid instability has the 
property, that it can generate an almost unlimited amount of resistive friction between 
the electrons and the waves. 

Now, let us imagine the two plasmas with opposite magnetic fields, Bq, approach 
each other. The pressure, between them, p = Bq/Stt, is dissipated at the rate Va/L 
by expansion due to fiow out the ends, and by force balance, must be replenished 
by compression due to dS/dt < 0. This compression normally continues until the 
Sweet-Parker thickness is reached. At this time, the plasma pressure in the layer 
is replenished by Ohmic heating rjj^ w r/i?Q/(47r(5)^ at the same rate at which it is 
depleted, {Va/L)p = (yA/L)(Bo/87r), by the adiabatic expansion. Therefore at this 
time the collapse dS/dt ceases. 

On the other hand, if the critical thickness, 5c-, is passed before the Sweet-Parker 
thickness is reached, the resistivity rapidly rises to generate an Ohmic heating large 
enough to balance the outflow adiabatic expansion at this larger distance (5c, and col- 
lapse ceases at this larger distance. 

For these conditions, the layer thickness is known, and Vr is determined by the 
mass conservation equation (2) alone, 



Thus, reconnection can become much faster than the Sweet-Parker rate based on 
Spitzer resistivity. Under solar flare conditions, Kulsrud (1998), it can become as 
much as a factor of a thousand faster. The resulting reconnection time can be reduced 
to a few hours, perhaps an order of magnitude longer than the observed energy release 
time in solar flares. 



5 Petschek Reconnection with Anomalous Re- 
sistivity 



In the second section it was shown that for a constant resistivity, Petschek's L' pa- 
rameter must be equal to L, so that Petschek's reconnection rate reduces to that of 
Sweet-Parker. However, if r/ is anomalous, enhanced by wave interactions, it can be 
very sensitive to the current density. The original problem with Petschek reconnection 
was that the external fleld dX x = L' was only slightly smaller, by a factor of 1 — L'^/L^, 
than its strength at a; = 0. But even this slight change in the resulting current density 
can lead to a finite and even large change in the resistivity rj. Taking this into account, 
one finds that equation (11) becomes 



Vr = 4Va. 



(18) 



dB 
'dt 



Att5L' 



Voc 



rj'c 




{rjo - r)')c 
Att5L' 



(19) 
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in which we have neglected L'"^ jl? and any shght difference between b and b' . r]o is the 
resistivity at x = 0, and rj' that at x = L'. Solving for By = {d9/dt)BQ as before, with 
this different value of dO/dt, and using it in Petschek's formula for the shock velocity, 
we find that with variable resistivity, 

V = = 1 Bp {vo-v')c 

^ Vkv^ 47r(5 ■ ^ ' 

Taking S = 6c, we have 

Now, to estimate the value of this revised reconnection velocity, we assume that rj 
is linear in j for j > jc, so that 

df] 1 Bo B' df] 1 L'^Bodf] 



Combining this with the mass conservation relation for the L' layer, 

Va L'' 



(23) 



(24) 



we obtain 

Yl = C d7] 

Fj VaL^ (47r)2 dj ■ 

This result can be written in a more familiar way by assigning a maximum value, 
T]* to T] and assuming that rj = rjspUzer for j < jc, and r] = rj* at j = 2jc. Thus, 

? = ^- (25) 
Prom jc = Bq/ AttSc equation (24) reduces to 

Vn fSclV^\ 



(27) 



Va \LS* 

where S* is the modified Lundqvist number based on rj* 

VaL 
77*c/47r 



Numerically, rj* comes from an electron wave collision rate equal to the electron 
plasma frequency uipe- Under typical solar flare conditions, Kulsrud (1998), rf* 

zer, SkTud 

10"^ (28) 



Va 

One can carry out similar estimates for the magnetosphere-solar wind interface and 
one finds from equation (26) that 

^ « 10° (29) 

Va 
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6 Conclusions 



We have shown or stated that: 

(1) . In general reconnection situations, L and are determined globally, while d 
and Vr are determined locally. 

(2) . For constant resistivity, the length of Petschek's diffusive layer is not a free 
parameter, but is determined by the condition that By be regenerated at the same rate 
as it is being dissipated by down stream flow. 

(3) . Constant resistivity gives L' = L, which makes Petschek's reconnection rate 
equal to that of Sweet and Parker. 

(4) . If the Sweet-Parker thickness S = L/s/S is thinner than the critical thickness 
Sc at which anomalous resistivity sets in, then the Sweet Parker reconnection rate 
becomes 

Vr = ^Va, (30) 



a rate that can be very much faster than their reconnection rate based on Spitzer 
resistivity. 

(5). In the case of anomalous resistivity the regeneration rate of By in Petschek's 
theory is much larger, and the Petschek's rate becomes faster even than the Sweet- 
Parker rate with enhanced resistivity. It is given by 

Vn ^*oM'/3 ^^^^ 



Va \LS 

where S* = LV^/(77*c/47r) is the Lundqvist number based on the maximum possible 
resistivity r]*. Note that it has a cube root dependence on this maximum resistivity, 
rather than a logarithmic dependence on the Spitzer resistivity which is the often 
quoted expression for Petschek reconnection. In spite of this, in many cases there is 
not a large numerical difference in the two results. Formula (31) gives an equally fast 
reconnection rate, and is more in tune with the true physical processes. 

(6) . A test for whether the anomalous resistivity rate equation (18) or (31), rather 
than the classical Sweet-Parker rate, equation (7), is applicable is: First, compute the 
Sweet-Parker thickness 5s p, of the reconnection layer Ssp = L/y/S, and compare it 
with the critical thickness Sc = Bo/{4TrneVi/c). If Ssp < Sc, then use the anomalous 
equation (31) for Petschek reconnection, or the anomalous Sweet-Parker equation (18), 
whichever is faster. 

(7) . In nearly all cases on the galactic scale, Ssp is larger than or at least com- 
parable to Sc, so the Sweet-Parker result gives the correct order of magnitude for the 
reconnection rate. This is almost always too slow to be of interest, so one concludes 
that reconnection on the galactic scale is hardly ever really important. 
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L' 

Figure 5: Sketch of the Petschek diffusion layer for the By field derivation. At x — L' , 
zero nearly up to y — 5' . 5{x) is a line of force. 

Appendix 

In this appendix we justify the intuitively described equation (11) for the evolution 
By, by a more precise derivation. In the intuitive derivation, it was assumed that the 
lines flowed into the reconnection only by resistive merging, and the eff^ect of plasma 
flow on them was ignored. Also, the merging velocity was taken proportional to the 
external field, through its effect on the current density. Further, the thickness of the 
reconnection layer was assumed constant between x = and x = L'. 

Consider an x dependent thickness S{x), as shown in figure 5, taken large enough 
that at y = 6{x),Bx ~ B{x) = Bo{l — x'^/L'^), the value for the external field. Also 
take 6{x) to follow a line of force. We have Ohm's law in the layer, 

- cE, = -VyB^ + V^By - r]cj, . (Al) 

In a steady state Ez is a constant. 

Integrate —cEz at a; = 0, from y = to y = Sq = 6{0). See figure 5. Since By = 
at a; = 0, we have 

fSo fSo /-So 

c 



/ -Ezdy = / -VyB^dy + rjc -jzdy, (A2) 
Jo Jo Jo 

or 5 5 5 



or 



roQ fon ^„ 

/ -E,dy= / -VyB^dy + ^Bo. (A4) 
JO Jo 47r 



Note that —Ez,—Vy, and —j^ are positive. 

Correspondingly, integrate —E^ at a; = L', from y = to y = S' = S{L') 



c f -E,dy = f -VyB^dy + f v^Bydy + ^Bo(l - L'^/L^). (A5) 
Jo Jo Jo 47r 
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At large enough S, we have 



-Vy{0,So) = VR, (A6) 

and 

since y = 5{x) is in the infinitely conducting region, and 

- cE, = -Vy{0, 5o)Bx{0, So) = -Vy{L', 6')Bx{L', S'). (A8) 

But, because is constant for all x and y, the left hand side of equation (A4) and 
(A5) arc equal to —cEzSq and —cEz6', respectively. Therefore, subtracting the right 
hand side of equation (A5) from S'/Sq times that of equation (A4), we get 



5' 

So -JO 



/ -VyB^dy- -VyBxdy+-^ —Bq-Bo(1-—) = -VyB^dy. (A9) 
JO Jo 47r Oq Jo 



The third term on the left hand side, the r] term, is the non uniform merging term 
quoted in the text. In fact, we can write it as a{ric/ A'k){Bq / doj^V^ / L?) where a is a 
constant of order unity and equal to 1 + {l/2)d'^6/dx^, which includes the change of 
S{x) with respect to x. The right hand side represents the down sweeping term which 
destroys By. 

We now argue that the sum of the first two terms on the left hand side is negative, 
so that the expression in the text for By should be an inequality rather than an equality. 
Since this is partially compensated by the fact that the regeneration term is slightly 
larger, by the factor of a > 1, we work with equality in the main text. 

First, the mean value of — in the first term, weighted by B^ is about 2/3 times 
—Vy{5o) = Vr because both —Vy and B^ are linear in y. Thus, 



S' r^o 6' 2 

-/ -VyB^dy^ --Vr^; (AlO) 

do JO Oq 6 



where 



^= t° B^dy, (All) 

JO 

is the Bx flux up to 5q. On the other hand, we see from figure 5 that the entire Bx 
flux at y = L' is near 5', so here the Bx averaged value of Vy is Vr/{1 — L'^/L^) k, Vr. 
Thus the sum of the first two terms in equation (A9) is essentially, 

where is the Bx flux at y = L' . Further, we also see, from Figure 5, that ^/^' is 
greater than il) because ij)' contains the flux (because 5{x) is a line of force) plus the 
flux through y = 0, between x = and x = V . Therefore, if (5'/(5 < 3/2, a plausible 
assumption, since L' <C -L, then the sum in equation (A12) is negative. In addition, 
the factors V'VV' > 1 ^'^d 1/(1 — L'^/L^) > 1 reinforce the conclusion that the sum is 
negative. 
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Now, dropping the first two terms in equation (A9) whose sum has been shown to 
be negative, we obtain the inequahty, 

a^^o^ > J^' -v^Bydy « Va6' < By > (A13) 

Comparing this with equation (11) of the text we see that, if a were one, equation 
(15) would yield an overestimate of By. Since a is certainly of order one, it is safe to 
take equation (15) as an estimate of By. 

Note, that equation (A9) is approximately the rate of change of L'SBy since 

£ dx dy^ = J J dxdy^ = d{x)E,\^' , (A14) 

and the latter expression can easily be identified with the difference of the two integrals 
given in equations (A4) and (A5) with the extra S'/Sq factor. 
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